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We consider searching for light sterile fermions and new forces by using long baseline oscillations of 
neutrinos and antineutrinos. A new light sterile state and/or a new force can lead to apparent CPT 
violation in muon neutrino and antineutrino oscillations. As an example, we present an economical 
model of neutrino masses containing a sterile neutrino. The potential from the Standard Model 
weak neutral current gives rise to a difference between the disappearance probabilities of neutrinos 
and antineutrinos, when mixing with a light sterile neutrino is considered. The addition of a B — L 
interaction adds coherently to the neutrino current potential and increases the difference between 
neutrino and antineutrino disappearance. We find that this model can improve the fit to the results 
of MINOS for both neutrinos and antineutrinos, without any CPT violation, and that the regions 
of parameter space which improve the fit are within experimental constraints. 

o 

q I. INTRODUCTION 

Non-collider experiments and astronomical observations have given us our first hints of physics beyond the Standard 
C$ Model, via the discoveries of neutrino oscillations, dark energy, and dark matter. The implication of these discoveries 
for fundamental physics is still unknown. The energy scale of the new degrees of freedom giving rise to neutrino 
oscillations could be as high as 10 16 GeV, as in Grand Unified theories, or as low as 0.05 eV, as in Dirac neutrino 
mass models. Even more mysterious is the nature of dark energy and dark matter, and the associated energy scale 
^ H or scales. If the new physics is light, it must be very weakly coupled to the Standard Model, or it would already 
have been discovered. Neutrino oscillation measurements offer an unmatched portal into any new nonstandard sectors 
containing light fermions, because neutrinos can mix with neutral spin 1/2 particles, and because oscillations over 
O-i long baselines are extraordinarily sensitive to extremely tiny effects. 

The long baseline experiment MINOS, which uses similar near and far detectors to reduce systematic errors, 
qj has observed the disappearance of both muon neutrinos and muon antineutrinos in the far detector PJHJ. The 
antineutrino data comes from a 7% antineutrino contamination of the beam and is severely statistics limited [5]. 
I— ~~ 1 Nonetheless, it is interesting to note that the antineutrino disappearance rate is larger than the rate expected from 
£^ neutrino disappearance by almost 2 sigma [4j. Currently MINOS is running in antineutrino mode, offering a unique 
opportunity to precisely measure the parameters governing the long baseline oscillations of muon antineutrinos. It is 
£SJ the purpose of this paper to offer a framework for searching for new physics in the antineutrino data. Most recent 
analyses [U [5HE1 15HTU] of anomalies in antineutrino data have focussed on CPT violation in the neutrino mass matrix 
[TTHl6| . However there is no theoretical motivation for CPT violation, and CPT violation requires Lorentz violation 
which is complicated to incorporate in a complete theory that is consistent with other data. In contrast, in this paper 
P\| we propose a simple, renormalizable, Lorentz invariant field theory, which is consistent with other experiments and 
which allows a significant difference between muon neutrino and antineutrino disappearance in the MINOS experiment. 
Many other papers have analyzed the consistency of neutrino oscillation data with sterile neutrinos [T7Tf2"2"] and 
' . I new forces |23j . Our model differs from those considered previously in three ways. First, we consider a relatively long 
range but weakly coupled new force for which the size of the matter effect can be considerably larger than the usual 
MSW effect, in a theory which is nevertheless consistent with precision electroweak constraints. Second, we consider 
smaller neutrino mass squared differences that primarily affect long baseline experiments for muon neutrinos. Thirdly, 
many previous analyses were concerned with the effects of sterile neutrinos on electron neutrino oscillations, while 
we are not attempting to address any anomalies involving electron neutrinos and are considering a model where the 
sterile neutrino has no electron neutrino component. 
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II. REVIEW OF THE STANDARD PICTURE OF NEUTRINO OSCILLATIONS 

A standard picture of 3 flavor neutrino oscillations has been successful in explaining phenomena observed by many 
diverse long baseline experiments [24, 25\. The e, /i, and r flavor eigenstates are related to the mass eigenstates via 
a unitary transformation, parameterized by three angles, neglecting a CP-violating phase: 

V e \ (l \ / C 13 S 13 \ / C12 Sis \ / VI \ ( vx \ 

V» = C 23 S23 1 -S 12 Ci2 \\ V 2 = U LNS \ »2 ] , (1) 
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where c^- = cos , = sin fty . A good fit to long baseline neutrino oscillation data may be found [26 for angles 

sin 2 20i2 = 0.87 ±0.03, 
sin 2 2<9 13 < 0.19 (90% C.L.), 
sin 2 20 23 = l.Otg.,,8, 

and neutrino mass squared differences: 

Am 2 2 = (7.6 ± 0.2) • 10~ 5 eV 2 /c 4 , 
|Am 2 3 | = (2.4 ±0.1) • 10~ 3 eV 2 /c 4 . 

Note that in this picture, for L/E < 1000 km/GeV, the smaller mass squared difference Am 2 2 gives oscillation proba- 
bilities which are always less than a percent, and which may be neglected compared with experimental uncertainties. 
Therefore for L/E < 1000 km/GeV, the observed oscillations are only sensitive to Am| 3 and, due to the small size of 
013, are primarily muon flavor into tau flavor. 

In matter, the oscillation parameters of neutrinos are modified due to the weak potential, a phenomena known as 
the MSW effect 27-29 . In the standard picture, only oscillations involving v e from a charged current weak potential 
are affected, as the and v r neutrinos only experience a flavor universal, neutral current weak potential in matter. 

The MINOS experiment studies the oscillations of both muon neutrinos and muon antineutrinos at a distance of 
735 km, over an energy range from 1 to 50 GeV. The survival probability as a function of energy for muon neutrinos, 
P{ v n v f/)i is well parameterized by Eqs. [2] and [3j The survival probability for muon antineutrinos, Pip^ — > is 
only marginally consistent with the standard picture, and is not well characterized by the oscillation parameters in 
Eqs.[2]and[3j The marginal agreement of the neutrino and antineutrino oscillation data at MINOS has raised interest, 
particularly because in vacuum, if CPT is a good symmetry, then the two oscillation probabilities are equal: 

CPT => P(i/ M -> = P(P M -> , (4) 

and in the standard oscillation picture, the presence of matter does not significantly change this relation at the MINOS 
baseline. 
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III. STERILE NEUTRINOS, A NEW FORCE, AND MINOS 

The MINOS antineutrino disappearance data has led several authors to investigate the possibility that CPT vio- 
lation may have been observed at MINOS. However, because the neutrinos are passing through matter, alternative 
explanations include a nonstandard interaction between the neutrinos and the matter which distinguishes between 
neutrinos and antineutrinos, a new kind of neutrino with a different weak charge, or both. 

In this work we discuss a model that can generate a difference between muon neutrino and muon antineutrino 
disappearance probabilities in matter at the MINOS energies and baseline. The basic framework is simple. We 
extend the Standard Model to allow the active neutrinos to mix with a sterile neutrino that has no Standard Model 
weak interaction. Because a sterile neutrino has a different weak charge than the active neutrinos, in matter muon 
neutrinos and muon antineutrinos will have different oscillation parameters, although this effect is small and only 
marginally detectable at the MINOS baseline. 

In Ref. [30], two of us proposed a new U (1) gauge force, coupled to the difference between baryon and lepton numbers 
(B — L), which can enhance this matter effect, in order to account for the anomalous results of the MiniBooNE and 
LSND experiments. We showed that provided the new boson was very weakly coupled (g^ 10 -5 ), such a boson of 
mass greater than w 0.2 keV could escape all laboratory constraints. If such a weakly coupled boson is lighter than 
a few MeV, it can provide a potential in the earth's crust which is much larger than the standard MSW potential. 
Furthermore, in Ref. [31j we showed that a light vector boson could have a chameleon nature, with very different 
effective mass in extreme environments, allowing it to escape astrophysical and cosmological constraints. Due to the 
very different energies and baselines of LSND, MiniBooNE, and MINOS, a minimal model which gives an observable 
anomalous neutrino/antineutrino oscillation difference at MINOS would not necessarily significantly affect LSND and 
MiniBooNE. 

The best fit points that we will find for the model from the MINOS neutrino and antineutrino oscillation data have 
a contribution to the potential from the B — L interaction which is approximately 4 orders of magnitude smaller than 
the one required to reconcile the shorter baseline MiniBooNE and LSND experiments. This allows the B — L coupling 
g to be lower by 2 orders of magnitude for the same B — L vector boson mass. Consequently, the constraints on the 
B — L force discussed in Ref. |30| do not limit the model in describing neutrino/antineutrino oscillation anomalies at 
MINOS, and in general at long baseline oscillation experiments. 

As in the usual seesaw model of neutrino masses [32H41] . we consider two types of neutrinos - neutrinos i>i with 
lepton number L = 1 which carry Standard Model weak interactions and neutrinos JVj with L = — 1 which are sterile 
under the Standard Model. In general the number of sterile neutrinos is not fixed by terrestrial experiment. We 
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allow for both lepton number conserving and violating mass terms, with both types of mass terms required to give 
active-sterile mixing. The allowed lepton number conserving mass terms, which arise from a Yukawa coupling to the 
Standard Model Higgs, are 



rriijViNj + h.c. . 

We also have lepton number violating Majorana masses for the sterile neutrinos 



MijNiNj 



+ h.c. 



(5) 



(6) 



The Majorana terms are gauge invariant under the Standard Model. With gauged B — L, they may arise from an 
expectation value of a B — L charged scalar. Unlike in the usual seesaw model, we assume that at least one of the 
sterile neutrinos is very light. For simplicity of analysis we will assume any other sterile neutrinos, such as would be 
required to cancel a gauge B — L anomaly, are heavier and more weakly mixed and will integrate them out. We thus 
consider an oscillation picture involving 4 light neutrinos- three active and one sterile. 

Neutrinos are produced and detected as one of the three active neutrino flavors, which propagate as a superposition 
of all 4 mass eigenstates. A 4-by-4 unitary matrix U transforms between the flavor and mass eigenstate bases. In 
general, four flavor neutrino oscillations are affected by 6 mixing angles and three CP violating phases. We make 
some assumptions to reduce the number of free parameters and simplify the analysis. We assume the neutrino masses 
to be hierarchical, with the heaviest mass eigenstate mostly sterile. We assume the sterile neutrino is a mixture of 
the 2 heaviest states only. We can then write the mixing matrix as 



U = U a U s 



(7) 



where U a is a block diagonal matrix 



U a 



Ulns 



(8) 



representing mixing between the three active neutrinos and 
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describes mixing between active and sterile neutrinos. In the flavor basis, the Hamiltonian is 

M m is the mass matrix in the mass basis and Vf is the potential matrix in the flavor basis, 

V f = diag(t/ cc - V nc - V B _ L , -V nc - V B _ L , -V nc - V B _ L , V B _ L ). (11) 

The potential V cc arises from the weak charged current, V nc from the weak neutral current, and V B -l from a new 
B — L gauge interaction. The values are: 

V cc = ±V2G F N ei V nc = ±%N n , and V B - L = ^N n , (12) 

where N e is the electron density, N n is the neutron density, g is the B — L coupling constant, and my is the mass 
of the B — L vector boson. With a reasonable matter density of p rts 2.7 g/cm 3 , the neutral current potential has 
magnitude 

V nc w 5.0 • 1(T 5 neV. (13) 

We will define the parameter V = \V nc + Vb-l 7 which will arise in the oscillation probability formula. The factor 
of 1/2 arises because the neutral current potential only affects active neutrinos. V in matter is positive for neutrinos 
and negative for antineutrinos. 

For simplicity, and because the MINOS range of L/E is not large enough to be sensitive to the smaller mass squared 
differences, we neglect the 2 smaller masses and take 



M m = diag(0,0, m,M). 



(14) 
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With this restricted mixing, U12 commutes with U s and M. m - Therefore, we can neglect Un in U a , since it does not 
contribute to oscillations. For simplicity, we also neglect the small angle 6*13. Then the electron neutrino does not 
significantly participate in oscillations at the MINOS baseline, meaning V cc is irrelevant and U a commutes with V/. 
We can further reduce parameters by assuming the mass m results entirely from mixing with the light sterile neutrino 
via the seesaw mechanism, receiving no significant contribution from the heavier neutrinos we have integrated out. 
The active-sterile mixing angle then satisfies 



tan (20 34 ) 



2^/]Wr 



M - m 

We can now write the Hamiltonian in an interaction basis where . 

H f = U a U' s H t U' s T Ul, 



with 



The diagonalized Hamiltonian is 



1 



U ' HiU ' T = 2E UsM ™ U > 



T 

m 1 ^ s 



V/- 



Hi = ^diag (0,0,m 2 ,A/ 2 
and the mixing matrix between the sterile neutrino and a linear combination of and v T is 

U' = 



1 














1 














cos 9 S 


sin 








— sin S 


cos 



In terms of the other parameters, 



where 



tan (20 s ) 



sin (26I34 ) 

COS (26*34 ) + 

AEV 



M 2 - m 2 ' 

Note that this mixing angle is larger in the antincutrino sector, where a is negative, and is maximal for 



cos (2 6*34 ) 



-AEV 



M 2 



\AEV\ = (M - mf 



(15) 
(16) 
(17) 

(18) 
(19) 

(20) 
(21) 
(22) 



This effect can give substantial conversion of muon antineutrinos into sterile fermions for a range of energies at the 
MINOS baseline. The eigenvalues of the Hamiltonian are 

, M 2 - m 2 



m 2 



M 



1 + a - y/l + 2a cos(2(9 34 ) + , 



M 2 



1 + a + + 2acos(26» 34 ) + a 2 



(23) 



The unitary matrix relating the interaction and flavor bases is 

1 

U a U' = 







c 23 s 23 c s s 23 s s 

-S 2 3 C 2 3C S C 2 3S S 
-S s C s 



(24) 



with c s = cos 6 S , s s = s'm9 s . Taking 6*23 — tt/A gives the simple mixing matrix 



U a U' s : 

The muon neutrino survival probability is then 

' rh 2 L 



/I 




Vo 











1 1 

V2 72 

1 1 

-i 



1 

V2 



\ 

s 
s 



(25) 



P(^ M ->• Vp) = 1 - cos 



, sm 



AE 



— sm 



, sm 



M 2 L 
AE 



— cos 



, sm 



, sm 



(M 2 



AE 



(26) 



The same formula applies to antineutrinos, although the values of m, M and 9 S are different because of the opposite 
sign of a in Eqs. [20] and [23] 
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IV. COMPARISON OF NEUTRINO OSCILLATION MODEL WITH MINOS DATA 

Armed with the muon neutrino survival probability, we can constrain the parameter space of the model using the 
results of the MINOS experiment. We will perform a combined fit of the model parameters m, V, and M using the 
MINOS i/fj, and P M survival data. We briefly discuss implications for the K2K [42] and SuperK [43] experiments, and 
the constraints from the CDHS results |44) . 

We perform a x 2 fit to the MINOS 



and data, simultaneously varying the three model parameters m, M, and 



V. The parameters are restricted to the ranges 



1(T 4 eV < m < 1 eV 
0.1 eV < M < 10 eV 
0.25 • 10" 4 neV < V < 10 neV 



(27) 
(28) 
(29) 



and each parameter is sampled from a logarithmic distribution (so that, for example, the same number of points are 
sampled for M in the ranges (0.1, 1) eV and (1, 10) eV). We use as our x 2 function 



X 2 = 2 iV oxp - iV obs + N obs In 



(30) 



which is more suitable for samples with small statistics [35]. To narrow down the range of parameter space, we use 
an approximation the x 2 value to estimate which regions of parameter space will provide good fits to the data. A 
sampling of approximately 2.5 • 10 5 points was sufficient to fill out the parameter space and identify the 68.3%, 95.5%, 
and 99.7% confidence level contours. A more detailed sampling was used to identify the best fit point. 

The MINOS data is divided into 24 total bins, 17 in the neutrino data and 7 in the antincutrino data. Fitting the 
three parameters m, M, and V yields a best fit point 



0.0394 eV, M = 0.157 eV, V = 2.01 ■ 10~ 4 neV. 



(31) 



The best fit point has a total x 2 = 24.8 = 1.24/dof. Note that the value of the potential at the best fit point is about 
8 times the neutral current contribution to the potential V, which is 0.25 • 10 neV. The mostly-active mass m at 
the best fit point is close to the corresponding value of (Am 2 ^) 1 / 2 in the standard picture, 0.049 eV. 
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FIG. 1: Regions within 68% (blue), 95.5% (red), and 99.7% (green) confidence limits, along with the best fit point 
(star), for pairs of parameters. The neutral current contribution to V, which is the lower limit, is the dashed line. 



In Fig. [TJ we show the regions within the 68.3%, 95.5%, and 99.7% confidence limits of the best fit point for pairs 
of parameters. In each plot, we have projected out the other parameter, taking the minimum x 2 over the parameter 
not displayed. The figures clearly indicate there are two regions of parameter space that are favored by the model. 
Both regions have m 2 near Am 2 tm for energies near the oscillation dip for both neutrinos and anti-neutrinos. 

The two regions of parameter space fitting the MINOS data have very different characteristics. One region has a 
nearly constant value of m, approximately the best fit value of m, and prefers larger values of M and smaller mixing 
with the sterile neutrino. This region contains the natural extrapolation of the model to the standard picture, where 
m 2 = Am 2 tm and M oo, which removes mixing between active and sterile neutrinos. 

The second favored region has larger values of V and smaller values of m. In the second region m and M are inversely 
related: as m increases in this region, M decreases. The allowed parameters produce an oscillation maximum for 
neutrinos at the MINOS baseline at about the same energy as the standard picture, while antincutrinos have a 
nonstandard looking pattern. In the second region the oscillation pattern for neutrinos can differ significantly from 
the standard picture in low density material and for neutrino oscillations at different baselines. 

Two other experiments with evidence for muon neutrino disappearance are K2K and SuperK. K2K uses a muon 
neutrino beam with baseline of 250 km and energy range from 0.5 GeV to 3 GeV. SuperK detects neutrinos produced 



by cosmic rays in the atmosphere and infers a baseline distance from their direction, with upward going neutrinos 
traveling much farther than downward going neutrinos. SuperK cannot distinguish neutrino events from antineutrinos, 
however most of their events are from neutrinos due to the larger neutrino cross section. K2K does not have any 
antineutrino data. SuperK data prefers oscillations of muon neutrinos into tau neutrinos over sterile states |45j . In 
our model muon neutrinos at SuperK should mostly oscillate into r neutrinos, since the mixing angle of neutrinos 
with sterile neutrinos is small at all energies and matter densities. However, due to the enhanced mixing angle of 
antineutrinos with sterile neutrinos in matter, for some parameters a significant fraction of the upward going muon 
antineutrinos at SuperK could convert into sterile states. 

We do not attempt an analysis of constraints from the K2K data, because the statistics are very limited, or of the 
SuperK data, because such an analysis requires access to details of the data. However we note that the first region 
of our fit looks very much like the standard picture at shorter baseline experiments such as K2K, where the effects 
of the potential are small, and in low density material, such as experienced by downward going neutrinos at SuperK, 
since in the first region the potential is small, mixing with the sterile state is small, and muon neutrino oscillations are 
dominated by m, which is near the standard atmospheric value. In the second region however the potential is larger 
and the neutrino oscillation pattern looks nonstandard at all distances other than the MINOS baseline. This region 
could conceivably be constrained by an analysis of SuperK data. For instance, the muon neutrino survival probability 
for downward going neutrinos at SuperK will be higher than the standard picture because matter effects are small 
in the atmosphere and because the vacuum mass squared difference between the neutrinos with large mixing angle 
is smaller. A SuperK analysis was done for mass varying neutrinos [46j |47] whose mass was different in the earth 
and in the atmosphere, and this analysis did produce some constraints |48j . In particular, a model with no muon 
neutrino oscillations in air was excluded at a 3.5 a level. Thus a similar analysis should be able to exclude the regions 
of parameter space with very small m. 

All of our fit points are consistent with short baseline tests of muon neutrino disappearance. The strongest such 
constraints come from CDHS [U] which is sensitive to mass squared differences larger than 0.3 eV 2 , and has the 
strongest mixing angle constraints for muon neutrino mixing for a mass squared difference of 2.5 eV 2 . For this model 
the we find the CDHS constraints are weaker than those of MINOS, due to the fact that the sterile mixing angle 
is quite small for larger mass squared differences. For comparison, we give the the effective mixing angle 9 c s for 2 
neutrino interpretations of muon neutrino disappearance at short baseline: 

sin 2 (2(9 0ff ) w sin 2 0„(1 + cos 2 6 S ) . (32) 

For the case where a is negligible, the fit points always have sin 2 (20 c ff) smaller than 0.13 for effective mass squared 
differences larger than 0.3 eV 2 . For effective mass squared differences larger than 2 eV 2 , sin (20 e ff) < 0.05, compared 
with the strongest CDHS constraint of sin (20 e ff) < 0.053. Nonnegligible a makes the effective neutrino mixing angle 
smaller. Thus the CDHS results do not provide any additional constraint. 

We also remove the B — L potential from the model, so that only the neutral current potential is present, and 
perform a two parameter fit to m and M. This fit finds a minimum % 2 = 28.1 = 1.34/dof, and lies within the 68% 
confidence level contour of the best fit point for the three parameter fit. The best fit values are 

V B -l = : m = 0.0420 eV, M = 0.309 eV. (33) 

and in good agreement with the MINOS neutrino and antineutrino data. 




FIG. 2: Muon flavor survival probability for the various models (solid red is the standard picture, dashed blue is the 
model with Vb-l, and dotted purple is the model without) against the data as a function of energy. 



In Fig. [2j we plot the predicted muon flavor survival probabilities as a function of energy at MINOS for the model 
at the best fit points with and without Vb-l against the data. For reference, we also include the survival probability 
for the standard picture. 

An analysis of the MINOS neutrino data has shown that oscillation into v T is preferred over v s . The model is 
consistent with this observation, as neutrinos oscillate primarily into v T . Antineutrino have larger oscillations into 
the sterile flavor than do neutrinos, and for the best fit parameters oscillations of muon antineutrinos into sterile 
neutrinos dominate oscillations in to v T at higher energies. In Fig. [3j we plot the r and sterile flavor conversion rates 
for our model at the best fit point for both neutrinos and antineutrinos at MINOS. In the same figure we also show 
the same data at the best fit point of the model with Vb-l set to zero. Note that the rate of conversion of active 
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FIG. 3: Predicted event rates with data as a function of neutrino energy, for the best fit point. The solid blue curve 
shows the \x flavor survival rate, while the dashed red curve and dotted purple curve show the rates for conversion 

into t and sterile flavors respectively. 



neutrinos into sterile is predicted to be small, due to the small vacuum active-sterile mixing angle #34. However the 
rate for active antineutrinos oscillations into sterile neutrinos is enhanced in matter. 

The MINOS experiment has also searched for oscillations into sterile neutrinos via measurements of the neutral 
current rate at the far detector [4"9"H5T] . They report a ratio R of observed to expected neutral current events. The 
effect of v e appearance, which can mimic a neutral current, is treated as an uncertainty in the analysis. The result is 



R = 0.99 ± 0.09(stat) ± 0.07(syst)(-0.08) 
R = 1.09 ± 0.12(stat) ± 0.10(syst)(-0.13) 
(34) 



E < 3 GeV 

3 < E < 120 GeV 



where v e appearance is assumed to be negligible, and the last number in parentheses shows the effect on R of assuming 
the maximally- allowed z/ e appearance. These results are obtained from a neutrino run, with only a 7% antineutrino 
component to the beam. We have also computed R for the best fit point, with and without a B — L potential, in our 
model. We have neglected the possibility of sterile neutrino scattering via the B — L gauge boson and mimicking a 
neutral current, because we are assuming the B — L coupling constant is small enough so that this is negligible. We 



find: 

E < 3, GeV, neutrinos, best fit Vb-l 
E < 4, GeV, antineutrinos, best fit Vb-l 

3 < E < 50, GeV, neutrinos, best fit Vb-l 

4 < E < 50, GeV, antineutrinos, best fit Vb-l 
E < 3, GeV, neutrinos, no Vb-l 
E < 4, GeV, antineutrinos, no Vb-l 

3 < < 50, GeV, neutrinos, no Vb-l 

4 < £ < 50, GeV, antineutrinos, no Vb-l- (35) 

The antineutrino data is split up at 4 GeV energy instead of 3 because the lowest bin for the antineutrinos goes up to 
4 GeV. Note that while the neutrino neutral current rate is in agreement with the reported MINOS data, there is an 
interesting possibility that the antineutrino run at MINOS could detect a significant depletion in the neutral current 
rate, particularly at higher energies. 
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V. CONCLUSIONS 



There is no compelling evidence for an antineutrino oscillation anomaly at the MINOS experiment. The data is 
statistics limited, and the current discrepancy is not statistically significant. However, this discrepancy offers the 
opportunity to discuss models that produce interesting experimental signatures at neutrino oscillation experiments. 
A neutrino/antineutrino oscillation anomaly can be explained by CPT violation, but it can also be explained by more 
conservative, well-motivated, and theoretically sound models. 

As we have shown in this work, the addition of a sterile neutrino can generate a difference in the neutrino and 
antineutrino oscillation probabilities due purely to the weak neutral current interactions, and the addition of a B — L 
gauge interaction will coherently add to the potential and enhance the oscillation probability difference. The model 
we have proposed can fit well the reported MINOS neutrino and antineutrino data, and matter effects are a nice 
framework to explore neutrino oscillations outside of the standard picture. 

As long baseline neutrino experiments enter an era of greatly increased baseline, intensity, and diversity, it is 
important to keep in mind that the history of neutrino experiment is full of surprises. Long baseline neutrino 
oscillation experiments using antineutrinos offer a unique and sensitive window to possible exotic sectors and forces, 
and it is important to keep a variety of possibilities in mind when analyzing the data. 
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